Current instability and diamagnetism in small-diameter carbon nanotubes 
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We investigate the electronic instabilities in carbon nanotubes of short radius, looking for the breaJi- 
down of the Luttinger liquid regime from the singular behavior of the charge stiffnesses at low ener- 
gies. We show that such a breakdown is realized in the undoped (3,3) nanotubes through the onset 
of phase separation into regions with opposite electronic current. The phenomenology derived from 
this regime is consistent with the formation of a pseudogap in the single-particle spectrum as well 
as with a divergent diamagnetic susceptibility, as observed in the experiments carried out in carbon 
nanotubes of small diameter. 
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Carbon nanotubes (CN) are among the best candi- 
dates for the development of devices in molecular elec- 
tronics. For this reason, there has been great interest 
during the last years in the investigation of their elec- 
tronic properties. As the electrons are constrained to 
move in the reduced dimensions of the tubule, CN con- 
stitute a paradigm of strongly correlated electron system 
[1-4]. Thus, it has been possible to observe signatures 
of Coulomb blockade in short nanotube samples, as well 
as signatures of ID transport characteristic of Luttinger 
liquid behavior [5,6]. 

Yet there has been also clear evidence of superconduct- 
ing (SC) correlations in CN attached to suitable contacts 
[7,8]. Supercurrents have been observed in the samples 
with SC electrodes reported in Ref. [7], providing evi- 
dence of the proximity effect in the CN [9]. Moreover, 
SC transitions have been measured in nanotube ropes 
attached to highly transparent contacts [10]. More re- 
cently, measurements carried out in CN inserted in a ze- 
olite matrix have supported the conclusion that strong 
SC correlations should exist in CN of very short diame- 
ter, leading to a transition temperature Tc ~ 15 K [11]. 

The experiment reported in Ref. [11] has provided evi- 
dence of the diamagnetic behavior of the CN at low tem- 
peratures. From inspection of the I-V characteristics, it 
has been also appreciated the appearance of a gap in the 
single-particle spectrum, which has been used to infer the 
value of Tc [11]. However, the lack of evidence for a true 
SC transition stresses the fact that the mentioned obser- 
vations refer to properties of the individual CN, which 
cannot be coupled by electron hopping as they only have 
a weak interaction with the zeolite walls [11]. In general, 
the tunneling amplitude between metallic nanotubes is 
what dictates the setting of 3D Cooper-pair coherence 
and the value of Tc [12,13]. The existence of a low-energy 
scale intrinsic to the individual CN demands therefore a 
reexamination of the role of the dominant correlations in 
nanotubes of short radius. 

The nanotube diameter d — 4. 2 ±0.2 A reported in Ref. 
[11] is closer to the value calculated for a (3,3) nanotube 
geometry, although the presence of (5,0) nanotubes in 



the zeolite matrix cannot be discarded [14]. It has been 
shown by using the local-density functional method that 
the (3,3) nanotubes have the same band structure of typ- 
ical armchair nanotubes near the Fermi level, with a pair 
of subbands crossing at two opposite momenta [14]. The 
case of the (5,0) nanotubes is different in that they have 
one more subband crossing the Fermi level, with angular 
momentum L — Q in the zigzag geometry. We will show 
that the (5,0) nanotubes have a extended low-energy 
phase with strong charge-density-wave (CDW) correla- 
tions, while the dominant low-energy instability for the 
(3,3) geometry is given by the phase separation into re- 
gions with opposite electronic current [15]. We will see 
that this is consistent with the existence of a low-energy 
transition inherent to the individual CN, allowing also to 
account for the phenomenology reported in Ref. [11]. 

Focusing first on the (3,3) nanotubes, the interaction 
processes between different low-energy branches crossing 
at Fermi points kp and —kp can be classified by attaching 
to them respective coupling constants gf'' [16] . The lower 
index discerns whether the interacting particles shift from 
one Fermi point to the other (i = 1), remain at differ- 
ent Fermi points (i = 2), or they interact near the same 
Fermi point (i — 4). The upper label follows the same 
rule to classify the different combinations of left- and 
right-movers, including the possibility of having Umk- 
lapp processes (j = 3) in the undoped system. 

We will adopt in what follows a bosonization approach 
to deal with the interactions in forward-scattering (FS) 
channels, while the effect of the rest of interactions will 
be considered through the analysis of their low-energy 
scaling. Thus, we first introduce an electron density op- 
erator Pri{k) for each of the linear branches of the sub- 
bands crossing the Fermi level, where the index r = L, R 
denotes the left- or right-moving character and the index 
i = 1,2 labels the Fermi point. It becomes possible to 
decouple the FS channels by passing to the combinations 

Pr±{k)^iy{pri{k)±pr2{k))/V2 (1) 

where u = — 1, -|-1 for r = L, R , respectively. The hamil- 
tonian for the FS interactions can be written in terms of 
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fields ^±(x), such that dx^±{x) = Pl±{x) + Pr±{x), and 
their conjugate momenta Il±{x): 

Hfs = 1 [dxJ2 {vjs{^s{x)f +VNs{da.^s{x)f) (2) 

The four independent velocities are given by 

vj± =vp + (I/tt) (fff ) ± gi'^ - {gi'^ ± ^)) (3) 

=vp + (I/tt) (^f ^ ± + {gi'^ ± fff ^)) (4) 

where vp is the Fermi velocity. In terms of these 
quantities, the renormalized velocities of the liquid are 
u± = y/VN±vj±, while the charge stiffnesses are K± = 

\/vj±/vN± [17]. 

In general, the significance of the interactions is given 
by their scaling behavior at low energies. The scaling 
equations of the gf couplings have been obtained in Ref. 
[16] for an armchair geometry. Here we improve nonper- 
turbatively the equations writing the scaling dimensions 
in terms of the K± parameters, and we introduce the 
pertinent modifications for the (5,0) geometry. The new 
equations read: 
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(3) Jl) , „(3)„(1) , ^(3)„(2) , J3)(2). 



-2grgr + gr'gr + grgr' + grgn (10) 

where / stands for minus the logarithm of the cniergy 
(temperature) scale measured in units of the high-energy 
scale Ec of the ID model (of the order of 0.1 eV). The 
couplings only arise in the (3,3) nanotubes without 
doping (which is the likely experimental condition when 
embedded in the insulating zeolite matrix) and the Up 
and ub couplings (to be defined later) apply only to the 
(5,0) nanotubes. The rest of equations not written here 
retain the expressions given in Ref. [16]. 



The scaling equations have to be solved with initial 
conditions accounting for the competition between the 
Coulomb interaction and the effective interaction arising 
from phonon-exchange [18]. The latter turns out to be 
attractive in the backscattering (BS) channels and repul- 
sive in the Umklapp channels [12,19]. Its effective cou- 
pling g has a strength » 0.3 — 0.9 times vp/n for a {n, n) 
nanotube [19]. This is comparable to the strength of 
the Coulomb interaction in BS and Umklapp processes, 
which we have taken as w 0.2 e^/n following Refs. [2] 
and [3]. Given that ~ 2.7vp, the Coulomb repulsion 
only becomes dominant in forward-scattering processes 
mediated by the potential V{q) ~ {2e^ / n) log(l + go/'?), 
where q is the momentum and n a suitable dielectric 
constant [20]. Under these initial conditions, the scaling 
equations develop an unstable flow below certain energy 
scale, at which the g'^^ and g'^^^ couplings enter a regime 
of large attraction while the g\ ' and Umklapp couplings 
are driven towards large repulsion in the (3,3) nanotubes. 

We have discerned the character of the electronic in- 
stability by looking for the susceptibility with the fastest 
growth as / — > oo, as well as for the possible breakdown 
of the Luttinger liquid parameters K± . We have checked 
that this latter instance takes place before any CDW or 
pairing response function starts to diverge in the (3,3) 
nanotubes. Their phase diagram is represented in Fig. 
1. The strength of the Coulomb potential V for the sam- 
ples reported in Ref. [11] has to be determined taking 
into account the screening effects from the large 3D ar- 
ray of CN inserted in the zeolite channels. These effects 
can be studied by generalizing the many-body approach 
devised in Ref. [21]. Thus, by performing a RPA anal- 
ysis incorporating the electrostatic coupling between all 
the CN in the array, we have found that the Coulomb 
potential ranges from V « 0.5 at 5 = 10 nm^^ to a 
saturation value V « 0.9 for q = 0.1 nm^^. Taking 

« 2.7vp and \g\ « 0.2-0.3 vp for the (3,3) nanotubes, 
we observe that the relevant regime corresponds to the 
phase Kj^ = 0, with very weak SC or CDW correlations. 

The vanishing of the Kj^ parameter corresponds to an 
instability in the balance between left- and right-moving 
electrons in the system. This follows from the fact that 
the vanishing of vj+ in (3) is what triggers the instabil- 
ity. A phase with = has been already described 
in the study of CN with purely repulsive interactions [4] . 
This instance has been also analyzed in Ref. [22]. In the 
case of CN of short radius, however, we have a nonnegli- 
gible strength of the phonon-exchange interaction. This 
leads to a different physical situation as the velocity vn+ 
remains finite, and the compressibility K^/u^ is nonvan- 
ishing at the point of the transition. We will see that the 
vanishing of vj+ marks actually the divergence of the 
n_|. correlations, opening a regime of phase separation 
into regions with opposite electronic current. 

The onset of the phase characterized by the vanishing 
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of Kj^ explains the apparent gap that has been measured 
from the I-V curves in Ref. [11]. The development of 
the gap shown in Fig. 2 of that Reference takes place 
through the change in the functional form of the I-V 
characteristics as the temperature is lowered, ressembling 
the behavior of a power-law dependence with increasing 
exponent. The same kind of behavior has to be found in 
the differential conductance dl/dV, which gives a mea- 
sure of the density of states near the Fermi level. 



one or two of the incoming modes belong to the subband 
with L = 0. These interactions may be classified with 
an additional set of couplings fj;''\ where the indices i 
and j keep the same meaning as for the gp'' couplings. 
Finally, there may be processes in which two particles 
with opposite L around the Fermi points 1 and 2 end up 
in the subband with L — 0, requiring the introduction of 
new couplings up and ub, for the respective cases with 
and without change of chirality of the particles. 
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FIG. 1. Phase diagram of (3,3) (left) and (5,0) (right) 
CN, in terms of the Coulomb potential V (in units of vf) 
and the elfective coupling of the phonon-exchange interaction 
G = 4:\g\/nvF- The thick lines are boundaries between the 
different phases characterized by the breaking of the param- 
eters K± and Ko at low temperature scales. In the (3,3) 
nanotubes, CDW correlations at zero momentum are dom- 
inant at the singular point, but with low strength repre- 
sented by the contours of constant CDW response function 

= 2, 4, 8, 16, 32, 64 (from top to bottom). For the (5,0) nan- 
otubes, three different regimes arise with strong correlations 
at the singular point of K- or Kq, given by the dominance of 
the 2fc_F CDW (white area), the 2fc_F spin-density- wave (light 
area), or the s-wave SC response function (dark area). 

In our description, the depletion of the density of states 
n(e) near the Fermi level can be evaluated with the usual 
bosonization methods, with the result that 



n{e) 



^{K+ + l/K++K.+l/K-~4)/8 



(11) 



The density of states follows a power-law behavior with 
an increasingly large exponent as vanishes in the low- 
temperature limit, as represented in Fig. 2. The plots in 
the figure show the development of the pseudogap at low 
temperatures. We observe that the shapes of the curves 
obtained from (11) are in very good agreement with the 
form of the I-V characteristics reported in Ref. [11]. 

Shifting now to the case of the (5,0) nanotubes, their 
behavior can be studied by making the parallel of the 
above analysis. One has to incorporate an additional 
type of density operators for the subband with angular 
momentum L ^ (and Fermi velocity v'p), which has 
its own renormalized velocity uq and charge stiffness Kq. 
Moreover, there are new interaction processes in which 
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FIG. 2. Plot of the density of states n{e) for V = 0.75ttvf 
and 4:\g\/nvF ~ 0.2 . The curves correspond to val- 
ues of minus the logarithm of the temperature scale 
I = 3.21, 3.225, 3.232, 3.235 and 3.236 . 

The processes given by the couplings up and ub have 
been overlooked in other analyses of the (5,0) nanotubes 
[23], while they are determinant in the enhancement of 
the SC correlations. Their scaling equations are given by 
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where A 



1 - 1/4:K+ - 1/4K_ - l/2Ko + a/f V^«F, 



vp/v'p and a = 2/(1 -\-v'p/vp). The dominant con- 
tributions in these channels come from phonon-exchange 
processes, so that the interactions up and ub enter a 
regime of large attraction at low-energies. The flow of 
the scaling equations gives rise now to a fast growth of 
the 2kF CDW, the 2fci? spin-density- wave, or the s-wave 
SC response function, leading to the phases shown in Fig. 
1. Recalling the value of the couplings for the experimen- 
tal samples, we see that the relevant region falls in this 
case into the phase with strong CDW correlations. 

An important point in connection with the phe- 
nomenology reported in Ref. [11] is that the (3,3) nan- 
otubes fall at low energies into a phase which is consistent 
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with the large diamagnetic signal observed in the exper- 
iments. This can be seen by recalling that the vector 
part A of the electromagnetic potential couples to the 
electronic current. Actually, around each Fermi point, 
the modes of the electron field can be arranged into a 
bispinor ^, whose hamiltonian is obtained from that of 
the low-energy excitations of graphene [24]: 

Hg = -ivF j (fr ^\t)<j-{V - i(e/c)A)*(r) (14) 

The ID projection to the low-energy excitations of the 
nanotube can be done consistently in the case of a trans- 
verse magnetic field, by chosing the vector potential as 
usual in the longitudinal direction of the nanotube [25] . 
After diagonalizing the quadratic form in (14), the ID 
hamiltonian becomes 

-ffiD ^""P j ^^'^^ ~ (e/c)^|| (<^)) 

{¥j,ik)^ji{k)-^lik)^Uk)) (15) 

where ip is the angle that supports the modulation of 
the longitudinal component A|| around the nanotube. It 
becomes clear that the vector potential couples to the 
charge asymmetry 11+ = pRi - pLi + Pr2 - PL2- 

When there is an enhanced susceptibility for the charge 
mismatch between left and right branches, a mechanism 
of screening of the magnetic field takes place, similar to 
the usual screening of the electric field from the total 
electron charge. The response function for the current 
n_|_, computed at the bosonization level, is 

iu,i.,m,i-.,-k)) = ^^^,^, (16) 

After integration of the density fields, we obtain a con- 
tribution to the free-energy density which, in the static 
limit, has the form -A\\{(f)A\\(Lp) vj,(e/c)'^ / K+u+. The 
susceptibility (16) leads therefore to a large diamagnetic 
response in the regime where 1/K+U+ diverges. 

The instability in the current is triggered by the van- 
ishing of the renormalized velocity wj+ at low energies. 
The diamagnetic response is therefore enhanced by the 
vanishing of and u+ at low temperatures. The be- 
havior of the susceptibility Vple/c)"^ /K+u+ as a function 
of the temperature is actually in qualitative agreement 
with the divergence of the magnetic susceptibility of CN 
reported in Ref. [11]. 

We end up with a picture consistent with the experi- 
mental observations reported in Ref. [11], assuming that, 
as supported by the measurements of the nanotube diam- 
eter, most part of the CN contained in the zeolite matrix 
have preferently a (3,3) geometry. We have seen that 
the transition reported in Ref. [11] may be interpreted as 
the breakdown of the Luttinger liquid regime, originated 
from the onset of phase separation into regions with op- 
posite electronic current. As well as for graphene [25], the 



existence of degenerate points in the spectrum is at the 
origin of the diamagnetic signal in the CN. The peculiar 
feature in the nanotubes is that their ID character leads 
to the splitting of the electron system into regions with 
alternating direction of the current. This implies neces- 
sarily the appearance of domain walls in the configuration 
of the current. The physical picture provides then a ID 
analogue of a well-known phenomenon in higher dimen- 
sions, by which the disordering effect of solitons spoils 
the propagation of the gauge fields. This is the role of 
the divergent fluctuations and consequent domain walls 
in the current, which lead to a phase transition in ID as 
that only requires the breakdown of a discrete symmetry, 
namely that of the current direction in the CN. 
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